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We evaluate all the electron-phonon couplings derived from the one-body and two-body 
electronic interactions, in both the adiabatic and extreme non-adiabatic limit, for a dimer 
with a non-degenerate orbital built from atomic wave functions of Gaussian shape. We 
find largely different values of the coupling parameters in the two limits, as well as 
different expressions of the corresponding terms in the Hamiltonian. Depending on the 
distance between the dimer ions, some of the two-body couplings are comparable, or 
even larger than the one-body ones. 

1. Introduction 

The models containing interactions of electrons and phonons have been the object 
of intense study in the recent years, mainly due to the discovery of high supercon- 
ductors {HTS) and colossal magnetoresistance manganites (CMR). When trying 
to apply the theoretical results to the real materials, however, one is confronted 
with the difficulty that a large degree of arbitrariness is present, stemming from 
two origins. First, different forms of interactions have been introduced, essentially 
on an ad hoc basis. Second, the value of the interaction parameters are estimated 
from the experimental data, which, however, can be interpreted only after a definite 
theoretical model, i.e. a set of electron-phonon interactions, has been assumed to 
be relevant for the material under measurement. Therefore, besides unavoidable 
quantitative approximations in fitting the data, also qualitative uncertainties are 
introduced a priori in the analysis. The purpose of the present work is to provide, 
in a simple but not trivial model, both rigorous prescriptions about the admissible 
electron-phonon interaction terms, and the quantitative evaluation of the interac- 
tion strength parameters. The model we consider is a dimer where the electrons 
occupy a non-degenerate orbital described by Wannier functions built from atomic 
orbitals of Gaussian shape. All the one- and two-body electronic parameters in that 
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model have been analytically evaluated in Ref.l. The electron-phonon interactions 
will be evaluated following the general definitions introduced in Refs.2 and 3. 

The structure of the paper is as follows: in Sec. 2 we describe the model and 
the general procedure of evaluation; Sees. 3 and 4 are devoted to the analytical 
evaluation of the coupling parameters in the non-adiabatic, and adiabatic limits, 
respectively, also defining the correct form of the corresponding terms in the Hamil- 
tonian; Sec. 5 discusses the numerical evaluation of the coupling parameters; Sec. 6 
is devoted to the conclusions. The Appendices contain details of the analytical 
calculations. 



2. The Model 



In a general electron-phonon Hamiltonian H — H^i + Hph + H^i-ph the interacting 
term Hei-ph originates from the development of the electronic part Hei to first 
order in the phonon-induced displacement of the ion positions. In our analysis we 
shall consider the terms due to the development of both the one-body andjtwo-body 
contributions to the electronic Hamiltonian, which, in standard notationQ, reads: 



'(nig- + n2g) + y^[t+X(ni^er+n-2-o-)](4aC2tT + jJ'.c.) + [/(ni|ni;+7j,2|n2i) 



-I- Vnin2 - 27^5*1 5*2 - JxyiStS^ + H.c.) + P{c{^c\^C2iC2^ + Ha 



The electrons are assumed to occupy a non-degenerate orbital described by Wannier 
functions built from atomic orbitals of Gaussias shape. 

2.1. Undisplaced Wannier functions 

The ions occupy the undeformed positions i?^ = (—a/2,0,0) and 1?° = (a/2,0,0). 
The unit vector ei2 ~ {R2 — -Ri)/a = (1,0,0) points from ion 1 towards ion 2; if 
not specified otherwise, e = 612. The position of the electron is r = {x,y,z). We 
associate to each site i = 1, 2 a Gaussian atomic-like orbital (f)i{r — Ri) . By defining 
N = (2/7r)^/'^r3/2, so that < (l),\(f), >= 1, they read: 

Mr - R°) = Ncxp j-r^ [(a; + a/2)^ + _^ z^] } , 

Mr - R") = Ncxp j-r^ [(.T - a/2)^ + + ^^j | . (2) 

Their overlap Sq = (0i|02) = exp(— r^a^/2) is non-vanishing. Then the Wannier 
functions ^'i, 5*2 can be written as: 

*i(r -R^r- Rl) ^ A(5o)0i(r - R\) + B{So)Mr - R2), 

*2(r- -R°,r- Rl) = B{So)Mr - Ri) + A{So)Mr - R^)- 
The requirement < >= 6ij yields the mixing coefficients A,B as: 



(3) 



1 



1 



B{So) ^ 2 



1 



1 



(4) 
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In the limit So 0, one has ^'i ^ (f)i {i = 1, 2), so that one can connect unambigu- 
ously each Wannier function to a well defined site. 

2.2. Displaced Wannier functions 

We shall consider only deformations = U\e\2 and U2 = ^2612 altering the length 
of the dimer, so that only their x-components are non vanishing. Their signs are 
defined, on both sites, with respect to the same unit vector 612, so as to have a single 
reference system. Symmetry requirements impose Ui = —U2 so that the relative 
displacement u = U2 — Ui = —2ui = 2u2- Even though we shall never introduce 
the exact time dependency in the following calculation, it is important to keep in 
mind that the displacements vary periodically in time with the frequency O of the 
phonon. The instantaneous positions of the ions are therefore defined as 



R^ 



(-a/2 + ?ii,0,0), ila = H2 + W2 = (a/2 + W2,0,0). (5) 



In the adiabatic limit, the orbitals on each site are assumed to adjust instantaneously 
to the displaced positions of the ions. Then, if we define 



!.i[r- = Afexp{-r2 [(x + a/2 - wif + j/^ ^ ^^j | ^ 



[x — a/2 — U2f' + : 



(6) 



^2[r- (il^ + Ma)] =^exp{-r2 
the Wannier functions 'fi (i = 1, 2)arc written as: 

Rl, Ml, U2) = A{ui, U2)Mr -R°i- ui) + B{ui,U2)Mr - -R2 - '"2), 

^2{Ri, R2, ui,U2) = B{ui, U2)Mr -Rl- ui) + A{ui,U2)Mr -R^- ^2). (7) 

In principle we have four coefficients A{ui,U2), B{ui,U2) to determine. We make the 
reasonable guess that, as for Ui the A's ( and the i?'s) have the same shape on 
both sites, then they keep the same shape also for Ui ^ 0. We shall define Ao,Bo, So 
aslimui,u2^oA{ui,U2) etc. As = {(I)2{u2)\(l)2{u2)) = 1 still holds, by 

defining S = exp{— (r^/2)[a — {u2 — wi)]^} the condition {^i{ui,U2)\^ j{ui,U2)) = 
Sij now yields: 



A{S) 



1 



B{S) 



1 



1 



(8) 



Namely, A{S) and B{S) depend on S{u) as Aq and Bq depend on So- 



2.3. The electronic parameters 

Given the Wannier functions, all the electronic parameters can be evaluated: this 
was done in Ref.l. To make clear our method of calculation, it is convenient to 
explicitate the one-body electronic interactions, corresponding to the local energy 
e and to the hopping amplitude t. 
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J ^*{r,R,,R2) 



2m 



+ Vi{r- Ri) + V2ir - R2) 



(9) 



t= / ^l{r,Ri,R2) 



_ + Vi{r- Ri) + V2{r - R2) 

Zm 



^2ir,Ri,R2)d^r. 



(10) 

In the formulas above, the potentials originating from the ion cores at the displaced 
positions Ri and R2 are: 



Vi EE V{r - Ri) = -e^Z 



V2 = V{r - R2) - -e^Z 



/ a 
(^+2 



a 



- Ul 



U2 



2 I 2 

y +z 



-1/2 



-1/2 



(11) 



where — e is the electron charge, and +Ze is the charge of the ion core. 

The local energy e (actually site-independent) can be decomposed into three 
terms, respectively corresponding to the contributions from the Laplacian kinetic 
operator (ey ) and from each one of the ionic potentials (ey^ , e^^ ) . The modulation 
of the term where both charge distribution and potential refer to the same ion gives 
rise to an Holstein-type coupling. The other term, where the charge around one ions 
feels the displaced potential of the other ion, we shall call "crystal- field" coupling. 
In the literaturei this latter term is usually assumed to be negligible but we shall 
see that actually is of the same order as the other ones, and can even become the 
dominant one. 

The modulation of the hopping.^ amplitude t gives rise to the so-called Su- 
Schrieffer-Heeger {SSH) interaction^. Similarly to e, it can be decomposed into 
kinetic (iy) and potential (t Vi,V2) contributions. 

We need to distinguish between the adiabatic (labelled ad) and the non-adiabatic 
(labelled no) limit in evaluating the electron-phonon interactions, because the in- 
tegrals have different kernels in the two cases. Indeed, in the adiabatic limit, when 

< t, the displacements affect both the potentials and the electronic Wannier 
functions, expressing the requisite that the electronic charge distribution adjusts 
itself instantaneously at the position of the ions. 

We shall schematize the opposite situation hH. > t, where the electrons are 
slower than the ions, as realized by the electronic charge distribution staying cen- 
tred around the undisplaced ion position, while the potentials are centred on the 
displaced ions. We shall call this the extreme anti-adiabatic limit. Not only the 
strength of the interactions, but also the form of the terms contributing to the 
Hamiltonian will turn out to differ in the two cases. 

The anti-adiabatic limit will be treated first, as the notation is easier to establish 
in that case. 
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3. Couplings in the non-adiabatic limit 

If the Wannier functions keep being centred around the undisplaced ionic posi- 
tions, neither ey nor t-^ does change, therefore no electron-phonon coupling orig- 
inates from them. The couplings derived from the two-body interactions are also 
identically vanishing in this limit, because they involve the Wannier functions and 
the inter-electronic Coulomb potential which are both insensitive to the displace- 
ments of the ions. 

The only non-vanishing types of electron-phonon non-adiabatic couplings are 
those arising from the variation of the potential contributions to e and t. 

3.1. The Holstein-type coupling go 

The Holstein-type coupling go is usually introduced as the site-independent am- 
plitude of a term in the Hamiltonian connecting the local charge with the local 
deformation: 

90^iniaUl+n2crU2)- (12) 



This interaction originates from the perturbation of the on-site atomic energy eo 
due to the ion motion in the non-adiabatic limit, namely: 



gl;^n,u, + Oiu^), (* = 1,2) 



(13) 



Our 5q*' is therefore different from the interaction defined by HolsteirH because 
he considered the interaction as due to the variation of the internuclear distance 
between the two components of each dimer attached at the nodes of a frozen chain. 
We use the "Holstein" label because this interaction has the same form as the 
traditional Holstein term. We also allow for a possible site-dependency of go \ to 
show later that this is not the case. 



For the ion at Ri one writes ^q^"* as: 

d 



(1) - I r "'^Vi 
gy = \ lim — — !- 

^" lui^o dui 



= I lim 



i{x,y, z, Rl, RlfVi{u^)dxdydz. (14) 



In the integral 



/oo 
■^i{x,y,z,Rl,RlfVi{u-^)dxdydz 
-OO 

-|).»,=1}V.[( 



a 

"+2 



ui ,y,z 



dxdydz, 
(15) 

let us change variable form x to p = x + a/2 — Ui, to shift the origin of x coordinate 
onto the displaced ion, without changing the shape of the Wannier functions. The 
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reason is to avoid the derivation with respect to Ui of Vi{ui) which is discontinuous 
in the integration range. The above expression changes into: 



{Ml [(P + Ui) ,y,z] + B(j)2 [{p-a + ui) , y, z]} Vi [(p, y, z)] dxdydz. (16) 



In this form the kernel has a continuous derivative with respect to Ui everywhere. 
Notice also that A and B still do not depend on Ui- To derive with respect to Ui 
Eq.^, we use = — Ui and notice that: 



lim ■ 



KiP + '^i^V^z) _ d(j)i{x,y,z) 



dx 



lim ■ 



yielding: 



d 2 
lim — {A(f>i [{p + u^) ,y,z]+ B4>2 [(p- a + u J , y, z]} 

u, ^0 U^ 



Therefore we can write: 



2*1 



dx 



dx 



2*1 



9*1 
dx 



2 / *i— J-y dxdydz. 
J-oo ox 



i{x,y,z) 
dx 
(17) 



(18) 



(19) 



By substituting *i and its derivative into Eq.p^ and by noticing that (^i | {x + a/2) 14 |0i) 
0, as the kernel is odd in a; + a/2, we obtain: 

/oo />oo 
{B^(t)l + 2AB(j)i(j)2)xVidxdydz + 2r^aB^ (j)lVidxdydz. (20) 
-oo J —oo 



The integrals in Eq.|20| are listed in Appendix 2. The final result is 

'a^r^ 



(1) 2rV2/7r 

?o = 

a 



B^ [Fo{2a^T^) - S^] + AABS 



Fn 



- S 



(21) 



That (7q is site-independent can be proved by considering that, for equal charges 
and displacement amplitudes, the energies = gQ^n.u, • e^j = 1, 2) on both 
sites must coincide, i.e. 



(1) (2) (2) 

'ui • ei2 = U2 • 621 -g^ 'u2 • ei2 

Now symmetry requires ui = —U2 from which gj^"^ — g^^"^ = go follows 



(22) 



3.2. The crystal-field coupling 

This term expresses the change in the energy of interaction between the charge on 
site i and the potential centred on the site j, that is: 



(*°|y,(HO + z.,)l*°> - (*°|1^,(J?°)|*°) ^ gi'jn^u, + 0(nf ) (23) 
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In the literature, there is some confusion on the form of this term, so we shall 
discuss it in some details, following the line of reasoning used to prove that gg^^ = 

(2) 

5q ' . Consider site 1, with charge ni. Its energy, after a displacement U2, changes 
by an amount E^^^ = g^^fUx U2 • 612- This can be considered as the quantity 
measured by an observer sitting on ion 1 and watching the ion 2 moved by U2 
. The equivalent measurement done by an observer on ion 2 watching the ion 1 
displaced by ui, yields E^'^^ = g^^^n2VLi • 621. Notice that each observer uses its 
own unit vector, pointing from the observer's ion towards the other, displaced ion. 
Formally, i?*^^^ can be obtained by performing a site label permutation on E^'^\ 
An external observer would find it more convenient to refer both quantities to the 
same reference system, ei2 say, yielding = —g^^^n-i ui •ei2. Assuming equal 
charge density and displacement amplitude in the two observations, one must have 
E^^^ = i?^^). It follows, dropping the charge densities: 

(1) (2) (1) (2) (1) (2) 

gl/u2 • ei2 = -flc/ ui • ei2 =^ gl/ui = -gl/u2 =^ gl/ = gy = g^f, (24) 

because the only displacements allowed by symmetry are such that ui = —U2- 
Therefore for the dimer as a whole one writes this term as 

9cf^i'f^i<TU2 + n2^ui). (25) 
Passing now to the explicit evaluation for site 1 we have: 

gi)^ = lim r m^{RlRl,x,y,zfV2{u2)dxdydz. (26) 

We change variables from x to p = x/a/2 — U2 and proceed in strict analogy to the 
evaluation of 5q^\ arriving at: 

^i-g^V2dxdydz. (27) 
By explicitating \E'i and its derivative we obtain: 

/oo 
(ylV? + + 2AB0102) xV2dxdydz 
-00 

/OO 
(ylV? - B'^(I>1) V2dxdydz. (28) 
-00 

^From the matrix elements evaluated in Appendix 2 we obtain: 



5fi = -2A (-\ J - [AFo{2a^r^) + 4BSFo{a^r^ /2) - ABS^ - AS^] . (29) 
■' \ a / V TT 
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3.3. The Su-Schrieffer-Heeger coupling 



This term is a non-diagonal inter-site coupling, due to the modulation of the hopping 
amplitude t. As the Wannier functions keep being centred on the static lattice 
positions, the difference in kinetic energy caused by the displacements is due only 
to the contribution to t from Vi and V2: 

J *i(r - i2?)[Fi(r -R°- u,) + V2{r -R^- iX2)]*2(r - 
- j ^i{r - R\)[Vi{r - R\) + F2(r - H^)]*2(r - R^)<fr 

= 712 ^(cLc2^ + CLC1^)(U2 - Ui) + 0[{U2 - uif]. (30) 
a 

Notice that, to preserve the invariance of the Hamiltonian under site permutation, 
the SSH coupling has to be odd under the same operation: 712 = —721 (see e.g 
Refs. 2,3). 

Let us distinguish the two contributions as ty^ and ty^ , defined as: 



t 



-f 

J — C 

-f 

J — c 



*i(u = 0)*2(" = 0) 
^{x + a/2 - ui)2 + y2 + 

«'i(u = 0)*2(w = 0) 



dxdydz, 



zdxdydz. 



(31) 



^{x-a/2-U2r+y' + z' 

To evaluate ty_^, let us define p = x + a/2 — Ui and use the symmetry relation 
Ui = —u/2, so that we can derive the kernel: 



dt 



Vi 



-f 

J — C 



9[*i(p, w, y, 2)*2(P. y, z)] 



du 

Next we notice that 

d[^i{p, u, y, z)4'2(p, u, y, z)] 



du 



Vi (p, y, z)dpdydz. 



lim 



du 



1 d[^i{p,u,y,z)'if2{P,u,y,z)] 

2 u^o dp 



so that finally we can write 

dt 

lim 



«— ►o du 



=-- r 

2 7-, 



dx 



■V^ dxdydz. 



(32) 

ia(*i*2) 

2 dx 
(33) 

(34) 



By performing similar manipulation on ty^{u2), but now with U2 = u/2, we arrive 
at: 

dt,, 1 r)f\T/.\T/_l 

(35) 



lim 



-[ 

2 7-, 



Summing the two contributions yields: 



'^''^''^Iv^dxdydz. 



dx 



712 



1 r 

"2y_< 



dxdydz 



a(*i*2) 



9a; 



^1 + :. 



dxdydz 



9(*1*2) 



9a; 



^2 = -A- + 3^. 

(36) 
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The first integral, A', after explicitating and ^2 becomes: 



1 /'OO 

x = - 

2 



dxdydz— [AB{(fi + 4>l) + [A^ + ^2)^^02] V^{x,y,z). (37) 



The derivatives 
9(0i(?!)2) 



4r^x</.i02, 



oa; ax 
when substituted into Eq.(|37|) yield 



-4r- 



dx 



-AT^ X 



(38) 



{Ai?[(0i|(x + |)i^i|^i) + (</>2|(x - |)yi|.^2)] + {A^ + b^){<p,\xvm} ■ 

(39) 

Similarly we get 



y = -2V^ \ AB 



y = — dxdydz 



9(*1*2) 



dx 



V2. 



.l|(x+ ^)F2|0l) + (<^2|(X - ^)V2|<^2)] + (^2 +i32)(0,|a,l/2|02)} . 

(40) 

The contributions from ((/)i|(a; + a/2)l/i|(/)i) in Eq.^ and (02|(a^ ~ a/2)V2|'/'2) in 
Eq.^ vanish due to the parity of the kernel. Substituting the matrix elements from 
Appendix 2 we get: 



7l2 



= { ^ [^0^ - (1 - 4a^r2)^^o(2a2r2)] | 



(41) 



Under site permutation a —* —a and A — )■ i? so that 712 = — 72iria as expected. 

In conclusion, in the non-adiabatic limit the electron-phonon Hamiltonian is 
given by: 



Trna 

el—phon 



.90 ^("-lo-Wl + ri2crU2) + ffc/ '^{niaU2 + n-2crWl) 

CJ i7 

+ 712^1(4.^2, +cLci,) ("2 -Ml). (42) 



4. Couplings in the adiabatic limit 



One can use the method introduced in the non-adiabatic limit to evaluate the 
couplings in the adiabatic case. However, one readily verifies that, in the explicit 
expression of the different electronic interactions in the adiabatically displaced state, 
u invariably enters in the combination a + u. Therefore the procedure of first 
deriving the integral kernels with respect to u, then taking the limit m — s- 0, and 
finally evaluating the integrals is equivalent to first evaluating the interactions for 



9 



M = 0, and then deriving them with respect to a. AU the parameters in Eq.|I] were 
exphcitly evaluated in Ref.l, so we shall simply derive them with respect to the 
dimer length a. 

4.1. The coupling term derived from e 

There is some confusion in the literature about the correct form of the electron- 
phonon Hamiltonian obtained in this limit from the variation of the local energy e, 
therefore we shall devote some space to discussing this point. 

In the adiabatic limit, Ui ^ in both the charge distributions and in the potentials. 
As the origin of the x-coordinate can be placed onto one of the displaced ions, 
this interaction couples the charge on site i to the position of site j through the 
modification of both the kinetic and the potential contributions. One has then 
to take into account the relative displacement of the ions. To obtain coupling 
terms from ev and ey it is essential that one uses the proper Wannier functions. 
If one instead adopted the local orbitals, such couplings would vanish. In the 
adiabatic limit, therefore, the overall e-derived electron-phonon coupling term in 
the Hamiltonian is: 

ffe^^ X! "'l'^(^2 - Ui) + g(^) ^ n2a{ui ~ U2) ■ (43) 
cr (7 

As (7e^^ = f/e^'' we can drop the site indexes, and write the total adiabatic contribu- 
tion from local energy terms to the electron-phonon Hamiltonian as: 

Kl-phon = 9e^{n2a - ni„){ui - U2) . (44) 

(7 

An expression similar to Eq.^ has been proposed in Ref.7. To convince the skeptical 
reader, in Appendix 4 we show, using the contribution from ty^ as an example, that, 
starting from the general definition of the coupling, one verifies Eq.^. 

4.2. Explicit expressions for the adiabatic couplings. 

4.2.1. Couplings originated from one-body electronic interactions. 

There are two type of couplings: one {g^), derived from the local energy: 

5,^g«-(-4^-t-5g (»,J = 1,2), (45) 

and one (712) from the hopping term (Su-Schrieffer-Heeger coupling): 

712 = 7i2v + 7i2V- (46) 

Below we list each contribution to the electronic interactions and the derived electron- 
phonon couplings. From the local energy on site 1 



10 



^V2 



-^e^ ( 2ry^ J [B^ + A'^Fo{2a''T^) + 2ABSFo{a''T^ /2)] , (47) 



one obtains the three contributions to 



9v 



2m 



(1 - ^2)2 



[2 (1 - a^T^ - S^)] , 



9vi 



-Ze^ 2r 



dA"^ 2ABS^ + B^S'^ 



du 



+ 



+ Fo(2a2r2) 



9u 



-Ze2|2r^/- (25) 



au a ^ ' 



5« = -Ze^ 2r 



dB"^ 2ABS^ + A^S^ 
+ 



du 
'2 



+ Fo{2a'T 



dA^ A^ 



du 



Ze' I 2r^/^ 1 (25) 



dAB AB 



du 



(l + a^r^)) 



(48) 



The two potential contributions depend on the site for which they are evaluated, 
but their sum, gy^ + , does not. 

For the hopping term, we can lump together the potential contributions which 
individually have no particular physical meaning. Then from 

t = t'^ + ty, 



tv = 



2m (1-52)' 



tv = -Ze^ {^^-^ [l + Fo(2a2r2)] +2{J^ ^ B'^)SFq{o?Y'' 12)} , (49) 

the Su-Schrieffer-Heeger coupling 7]^2 = 7i2V + 7i2V obtained as: 



7l2V 



2m 



aSV^ 



(1-52)2 



[2(l-52)-a2r2(l + 52)] , 



7i2y = -Ze^\ 4r 



dAB ^2 + ^2 
du a 



AB 
a 



52 + —s^ 



dAB AB 
du a 



2t^2 



Fo(2a2r 



11 



+ s 



(l + a^r^) 



9u 



Fo(a2rV2)J>. 



(50) 



Notice that, as the partial derivatives are linear in a., then both 7^2 v 7i2y 
change sign under site permutation, as expected from Refs.2,3. 

The electron-phonon Hamiltonian in the adiabatic limit has therefore the form: 



Hel-phon = 9e ^i^la - n2a){U2 - Wl) + 7l2 ^{cLc^^ + cLci<,)(W2 - Ui). (51) 



4.2.2. Couplings originated from two-body electronic interactions. 

Also the two-body electronic interactions U, V, J(= P),X of Eq. (1) give origin 
to electron-phonon couplings which, to the best of our knowledge, have never been 
considered in the literature up to now. The electron-phonon coupling originating 
from the modulation of the kinetic exchange in the t — J Hamiltonian has been con- 
sidered in Ref.l7. However, the kinetic exchange Jkin = 2t^/(C7 — V) is a composite 
quantity different from the direct exchange J^y = = P in Eq. (1). Here we list 
the values the two-body electron-phonon interactions have in our model. They were 
obtained by deriving with respect to a the two-body interactions evaluated explicitly 
in Ref.l. Notice that, as U{a) = J{a) + e^F/y^, then dU/da = dJ/da = dP/da 



dX 
da 



, r 



i-aVS) 



(l + 3g^) 
(1 - ^2)3 



1 + 2S'^ + Fo (a^r^) - 2(1 + S'^)Fo 



—e 



S/a 
(1 - 52)2 



a2r2 



S/a 



(1-^2)2 



2{1 + S^) 



Fn 



- 1 



dU 
da 



r 



-4ar252 



(1-52)3 

52/a 



52-Fo(a2r2)| 



(52) 



2 - 52 -h 25* -h 52i^o (a^r2) - 452^0 



2a^r 



2t^2 



(1-52)2 

■\S^ -Fq (a2r2) -h 4 



1 - 452 - (a^r2) 4i^o 



a2r2 



Fq 



a2r2 



V5 



da 



2c2 1 r 



r 

7^ 



AaT^S 
(1-52)3 

1/a 



3 - 52 - 85* - (7 - 552)Fo(a2r2) - 4(1 - 352)Fo 



(53) 

a^T 



(1-52)2 



2a2r252 



-l-452-FFo(a2r2) +AFo 



a2r2 
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25' 



5^- 2-^5' Foia'T') + 2S' 



Fo 



(54) 



dJ _ dP _ dU 
da da da 

All the above interactions change sign under site permutation. The terms they 
contribute to the electron-phonon Hamiltonian all have the same form, namely: 



dY 



- Uj) = 1, 2) 



(55) 



where Y = U,V, X, J and F{cl^, Cj^) = 1, 2) is the function of Fermi operators 
representing the two-body interaction whose amplitude is Y 



5. Results 

By definition, clS Tlio' IS the number of electrons per site (a dimensionless quantity), 
then the interaction parameters have dimensions [energy][length]~^ . To get them 
in energy units (eV) we shall measure the deformations Ui in units of the char- 
acteristic phonon length L — y/Ti/ (2r2M). We choose the phonon frequency such 
that hil — 0.1 eV, which is appropriate to HTS and CMR, and M equal to the 
mass of ^^O. Figs.l and 2 show the behaviour versus the dimer length a (in A) of 
the non-adiabatic coupling parameters 5o7 5c/,7i2 (in eV) evaluated by assuming 
two representative values F = 1.0, 2.0A~^ for the shape-controlling parameter of 
the Wannier functions. The most unexpected result concerns gcf. While usually 
neglected in the literatureB oja metallic systems, this coupling has been recognized 
as relevant to polar materialsB. We find indeed that, when F = 1.0A~^, gcf is larger 
than go for any a, and it becomes the largest parameter for a > 2.2A. For small a. 
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the SSH coupling is the largest. The strength of all the couplings decreases with 
a, less quickly for g^f, which, for large a, is still comparable to the values of the 
hopping amplitude t{a). 

When r = 2.0A~^(see Fig. 2) Qq is negligible for any a, 712 is large for small a but 
drops very rapidly to negligible values as a increases, while gcf keeps appreciable 
values for all a values. We can conclude that, in the non-adiabatic limit, the more 
localized are the orbitals, the more relevant is the role of gcf in relation to the other 
admissible couplings. 

Fig. 3 is the adiabatic counterpart of Fig.l. The Holstein type interaction is 
absent, being identically vanishing as discussed above. Here we find that g^ is 
always larger than the SSH interaction 712, and particularly for large a there is an 
order of magnitude difference between them. 

Fig. 4 shows the couplings derived from the two-body electronic interactions for 
the same parameters as Fig. 3. In general, their values arc smaller than those of 
g^ and 712, with the possible exception of dV/da. Indeed, that coupling arises 
from a physical mechanism not very different from the one originating i.e. the 
vibration of the charge on site j as felt by site i. Similarly to also dV/da decreases 
slowly with a, so that for large a those two are the only relevant couplings. It is 
worth stressing that, though U is larger than J, their derivatives coincide. Besides, 
the derivative of the interaction coupling X is a non-monotonic function of the 
lattice constant. Indeed, X exhibits a maximum for a ~ 1.8 and then sharply 
decreases to a negligible value. 

When r = 2.0A~^, as shown in Fig. 5, g^ and 712 have trends similar to the non- 
adiabatic case, with larger values for small a. As a increases, 712 quickly reduces to 
negligible values, while g^ decreases more slowly. Of the two-body coiiplings, only 
dV/da is non- negligible, and it has a strength very close (in absolute value) to (/<=. 
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Figure 5: Adiabatic coupling constants g^, ^12, dV/da (in eV) versus the dimer 
length (in A), for T 2.0 A^i. 



6. Conclusions 

We have presented the analytical evaluation of the electron-phonon coupling pa- 
rameters derived from both one- and two-body electronic interactions in a model of 
a dimer with a non-degenerate orbital built from atomic orbital of Gaussian shape. 
The approach wc have followed (proposed in general terms in Refs.2,3) consists in 
inserting the site displacements ui,U2 in the kernels of the integrals defining the 
one-body electronic parameters eo and t, and then considering the first-order terms 
in the expansion of 60(^1, U2]-,aJld piM-ii ^2)- In the adiabatic case this procedure is 
equivalent to the simpler onOT0fl^ll3 consisting in evaluating the electronic param- 
eters as functions of the lattice parameter (corresponding to the dimer length a) 
and then deriving with respect to a. In the anti-adiabatic case, however, the latter 
procedure can not be applied, because a and u enter independently the kernel of 
the integrals. From the physical point of view, the " a-derivative" method describes 
how the electronic parameters vary under a quasi-static (low frequency) change of 
the equilibrium position of the ions, as could be realized e.g. under pressure. In 
general the electron-phonon couplings are the effect of the ion oscillations at fre- 
quency f2, which causes the system to be alternatively compressed and elongated 
over a time « fi^^. It is this implicit dynamics of the displacements which allows 
for distinguishing between the adiabatic and non-adiabatic regimes, and forbids the 
use of the "a-derivative" method in the latter case. 

A novel result is the evaluation of the couplings originating from the two-body 
electronic interactions. We have shown that at least one of them, generated by the 
Coulomb repulsion between the charges on different sites, is comparable, or even 
larger, than the couplings derived from the one-body interactions. 

The quantitative results for the coupling parameters, even if agreeing in order 
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of magnitude with some estimates from experimental data (see, e.g. Ref.13) are 
obviously model-depending. However, their ratios should be more close to the real- 
ity. In particular, the obtained values of the various couplings, when compared to 
the values of the electronic interactions obtained from the same Wannicr functionsld 
suggests that, for dimer lengths comparable to the lattice parameters in HTS and 
CMR, only dU /da and dX/ da can be safely dropped, while neglecting any of the 
other electron-phonon interactions is a questionable approximation. 

Finally, we have been able to determine the correct form of the admissible 
electron-phonon coupling one-body terms in the interacting Hamiltonian. In the 
literature one sometimes finds proposals for such terms which are incompatible 
with our results. This is the case, for instance, of the electron-phonon Hamiltonian 
of Refs.14,15. 
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7. Appendix 1 



Here we list some results of frequent use in the calculations. We shall use the 
zero-displacement derivatives of A and B: 



OA 
du 

dB 
du 



M->0 



.(l + 5„)-3/2 + (i_5^)-3/2 



-(l + 5o)-3/2-(l-5o)- 



■3/2 



(56) 



(57) 



In the following we shall use the substitution t = x + a/2 — u\ in both ^\ and ^'2, 
yielding: 



^x{x,ui,U2,y,z) => '^i{t,u,y,z)=A{u)Ne-'^ I* +^ ^+B{u)Ne-'^ [{t-a+u) +y +z ]^ 

*2(x, ui,U2, y, z) ^ ^2{t, u, y, z) = B(w)iVe-r'[*'+'''+^'l+A(u)7Ve-r'[(*-»+")'+'''+^'l. 

(58) 

To write down their derivatives in the zero- displacement limit, it is convenient to 
revert to x = t — a/2 so that: 



d^i{x,y,z,u) 



du 

d^2{x,y,z,u) 



du 



dA 

= (j)i{x,y,z)— + (t)2{x,y,z) 



^9B ^ , 
9nx,y,z)-^ +(p2{x,y,z) 



du ^ 2> 



dA „„9, a, , 



(59) 
(60) 



We shall also use the other substitution t = x — a/2 — U2, yielding: 

*i(x,wi,W2,,y,2) ^ A(«)7Ve-^'[(*+»-")'+^'+^'l +B(«)iVe-r'[*'+^'+^'l, 

^2{x,Ui,U2,y,z) ^ B(u)Afe-^'[(*+'^-")'+^'+"'l + A(w)7Ve-r'[*'+^'+"'l. (61) 
In the limit of vanishing deformations it follows, after reverting to x = t -\- a/2: 



d^i{x,y,z,u) 



du 

d^2{x,y,z,u) 



<l>i{x,y,z) 



u->0 



du 

8. Appendix 2 



= (j)i{x,y,z) 



dA „^9, a, , 



.dB 

+ (/>2(x,y,z) — , (62) 



+ (j)2{x,y,z)^. (63) 



For convenience reasons wc shall list here various matrix elements, all evaluated 
in the limit of vanishing deformations, which enter the calculations. Those of 
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them which do not follow straightforwardly from standard properties of Gaussian 
integralslij are evaluated in details in the Appendix 3. By defining: 



we have: 

qp2 p4 2 c2 



2 '2(l-5g)' = -2lT^' 



2t^2\ 



(cl^m^,) - 2rV2/^, {(b,\V,\^^) = 2rV2/7r^^o(2a'r 



{(bi\xVi\(bi) = -arV2A, (02kV2|02> = aVy^, 
(01 1x^2 !</>!) = ^y|[-5o'+ (l-2a'r2)Fo(2a2r2)] , 
{^2\xVi\cl)2) ^~{cf>i\xV2\cl)i), 



aT 



So — Fq 



(011X^2102) = -(01^14102). 

Some other useful relations are: 

(0i|^^i|0i) - {HV2\h), {h\vM = {(bi\vM. 

We can directly evaluate 



(64) 

(«,J = 1,2) 

(65) 
(66) 

(67) 
(68) 

(69) 

(70) 
(71) 

(72) 



Indeed, due to the odd parity of the kernel, {4>i\{x + a/2)Vi|0i) — from which 
Eq.(^ follows. We can also prove that 



l\xV2\^2) = -ihlxVM, {cj)l\xV2\cj)l) - -(021X^1 102). 



(73) 



Indeed, if x changes sign, then Vi V2 and (pi — s- 02- Therefore the functions 
4>i{Vi + ^2)^2 and 01 V2 + 02^1 are even in x, so that: 



(0i|a;V2|02) + (02|a;14|0i) = (0i|a:(Vi + 1^2)102) = 0, 

/OO 
x{(t)lV2+(fiVi)dxdydz^Q, 
-OQ 

because the integrands are odd functions of x. 



(74) 
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9. Appendix 3 

Let us now evaluate the unknown integrals in Eq.(^0|). To do that, consider 
(^ilyil^-i) = / ^i{a/2,x,y,z)^Vi{a/2,x,y,z)dxdydz, 



and its derivative with respect to a/2: 



d{^i\Vi\-i'i 
d{a/2) 



9*1 



d{a/2) 
9*1 



Vidxdydz 



t2 dVi 



'^d{a/2) 



dxdydz 



Vi dxdydz — g^^ 



(1) 



(75) 



dia/2)_ 

where we use the Eq.(|T9|) and lim„j_^o 9Vi/97ii — ~dVi/dx — ~dVi/d{a/2). To 
evaluate the first integral we need also: 

Inside the integral, it is convenient to use d/dx, instead of d/d{a/2). As 
9*1 



= ~2T^ [A {x + a/2) 0i + B{x - a/2)02 



(77) 



then 



dx 



so that Eq.(75) becomes, recalling Eq.(|19| 
OA 



d{^i\Vi\^i 
d{a/2) 



= 2 



9(a/2) 



■^icjjiVidxdydz + 2 



dB 
d{a/2) 



i<j)2Vidxdydz 



/oo 
■i>i{x-a/2)(j)2Vidxdydz. 
-OO 

Developing = A(f)i H- B(f)2 and reordering, one arrives at: 



(79) 



{AB{^i\xVM + B^^2\xVM) = 
dB^ T 



9(*i|l/i|*i 



dA^ 



+ 



AaT^B 



2 o2 



9(a/2) 



(02|li|</'2) + 



AaT^AB - 2 



d{a/2) d{a/2) 
d{AB) 



(011^^1 1 01 



9(a/2) 



(80) 



This equation connects the unknown integrals on the left hand side to known ones. 
We need another relation, which is provided by similar manipulations on (*2|Vi|*2)- 



9(*2|V"l|*2) 

9(a/2) 



d 



= 2 



9(a/2) 

9*2 



*2(a/2, X, y, z)'^Vi (a/2, x, y, z)dxdydz 



* 



9(a/2) 



Vidxdydz 



*^ 



_dVi_ 

d{a/2) 



dxdydz. 



(81) 
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After substituting dVi/d{a/2) = dVi/dx the second integral is done by parts and 
one arrives at: 



d{'f2\Vi\'f2 



= 2 



d{a/2) 

^From the definition of ^'2 one has: 



d{a/2) dx 



Vidxdydz. 



(82) 



piyu no n/1 



dx 



= -2r2 (.T + a/2) + A(.T - a/2)02] 



(83) 



so that 



'"^ '"^ ' +^02+4r^A(.-a/2)^2. (84) 



9(a/2) 9a; 9(a/2)^' a(a/2) 

Substituting Eq.(p4[) into Eq.(^), developing ^2 = B(j)i +A(j)2 and reordering yields 
finally: 



2 a2 



dia/2) 



(02 1 14 1 02 



AaV^AB - 2 



dia/2) d{a/2) 
d{AB) 



d{a/2) 



1^1 102), (85) 



which is the second equation needed to evaluate the integrals entering the Eq.(|20|) 
for ^Q^^ For the explicit evaluation we need also: 

(^ilFil*!) 2rv/2/V [A^ + B'^FQ{2T'^a^) + 2ABSFo{T^a^ /2)\ , 

{^2\Vi\^2) = 2Vy/2l^ [B^ + A^FQ{2T'^a'^) + 2ABSFo{r^a^ /2)] . (86) 
It is convenient to subtract the Eq.(^o[) from Eq.(|85|) yielding: 

{A^ - B^) < (t>2\xVi\4>2 >= [< ^2\Vl\^2 >-< ^'ll'^^ll^'l >] 



d {A^ - B^ 
d{a/2) 



< 01 1^1 101 > + 



4T^a {A^ - B^) - 



d {A^ - B^ 
d{a/2) 



< 02 I Vl 102 >, 

(87) 



while their sum yields: 

ler^AB < 0l|xli|02 > +8T^ (A^+B^) < 021x^x102 >= 



d d(A^ + B^) 

'< *2|14|*2 > + < *i|14|*i >] - ^ < 0l|14|01 > + 



d{a/2) 



d{a/2) 
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.2 ,,2 , „2N diA'+B^) 



4T^a {A^ + B^) - 



d{a/2) 

d{AB) 



< 4'2\Vl\(j)2 > 



d{a/2) 



< <l)i\Vi\4<2 > . 



To proceed, one has 



< ^2\Vi\-^2 >-< ^-ilT^I^-i >= -r\/- (^^ - -B^) [1 - Fo (2r2a2)] 



< *2|Vl|*2 > + < *l|Vi|*i > = 



rJ-{{A^ + B^) [l + Fo(2rV)] +4ABSFo{T^a''/2)}. (89) 



The derivatives of interest are: 



dS 



5 (a/2) 



-2r^aS, 



(90) 



d {A^ + B^) _ 4T^aS^ d {A^ - B^) _ 2T^aS^ OAB _ T^aS (l + S^) 



a (a/2) (1-52) 



2 ' 



a (a/2) d{a/2) 

(91) 



2 ' 



d 



[Fo {2r'a')] =-[S'- Fo {2T'a')] , 



d 



a (a/2) 
Therefore we obtain: 



a (a/2) 



Fo 



S-Fo 
(92) 



2„2 \ "1 



d 



5 (a/2) 



[< *2|Vi|*2 > - < >] 



2T.j2f^ 
a (1-52)^/2 



{T^a'S^ [1 - Fo {2T^a^)] + (l - S^) [S^ - Fo {2T^a')] ] , (93) 



and 



d 



a (a/2) 



[< \E'2|Vi|*2 > + < ^ll^ll*! >] = 



a (1-52)' 



{-2T^a^S^ -{2-S) S"" (l - S^) - Fq {2r^a^) [l - S^l - 2T^a^)] 
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+ Fo 



(94) 



By substituting the equations above deduced and the other matrix elements into 
the Eq.(^7|) and Eq. (p8|) we obtain the unknown matrix elements as: 



{<l>l\xVi\<f>2) = 



2ar 



aT 



[54-Fo(2a^r2)], 



S-F: 



(95) 
(96) 



10. Appendix 4 



Let us evaluate the crystal field coupling for site 1 in the adiabatic limit, in order 

to verify Eq.Esl For the charge on site 1 the variation of the local energy defines 
(1) 

9cfad as: 

/OO rOO 
- OO J —OO 

^ 9i]lMu2 - ui) + 0[{ne - (97) 
Substituting t — x — a/2 — U2 and developing ^'i (u, t, y, z)^ to first order in u yields: 



0^1^ -2 



*i(0,t,y,z) 



d{u) 



V2it,y,z) • ei2. 



(98) 



By substituting t ^ x — a/2 — U2 we obtain: 



du 



u=0 



dA 
du 



u=0 



dB 
du 



-2T^A{t + a)Ne-^^^'+''^+y+''l 
Substituting *i and Eq.(p9|) into Eq.(|8|) yields: 



(99) 



{^l\V2\^l) + 



9(B2 



du 



(02 1 1^2 I 02) +2 



d{AB) 



du 



(01 1 1^2 1 02 



+ [Al{cj)i\ix + a/2)y2|0i> + 2^oBo(0ik^^2|02)] 



(100) 



The corresponding coupling for the other site is obtained by the site permutation 
1^2, 

^^c/L^ffgLE "2.(^1 -Z.2), (101) 

that is, with t = x a/2 — uii 



o^') -2 



^'2(0, i, 2/, z) 



5*2 



d{ux ~ U2) 



(102) 
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Performing the corresponding calculations we obtain: 

5c/L = - 4r4 AB{<j,^\x{V^ + ^2)102) + A\[{^^\xV2\^i) + {M^Vi\d,2))] 



+ B^{[{(j,2\xV2\(^2) + {4>i\xVM)]y (103) 

The difference between g^^j^^^ and g^^^^; vanishes because of the results reported in 
the Appendix 2, Eq.( |4|). 
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